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length = 64; //fast fourier transform length 



/ 



delta = .1; %input ( 1 Enter stop band tolerance: '); 
beta = 0.2; %input ( 1 Enter pass band tolerance: 
wp=8 ; 
ws=7 ; 

mu = 28; %CP length 

load channel_impulse_resp.dat; %load channel impulse respons 
N=19; % sirf filter length 

ch = channel_impulse_resp; 
mm=max (abs (ch) ) ; 
ch=2047*ch. /mm; 

NN=size (channel_impulse_resp) ; 
M=NN ( 1 ) ; 
for i=l:mu 

mv(i)=l; %initialize mv 

end 

C= conv (mv, abs (ch) ) ; 

[Y,I] = max (abs (C) ) ; %find max which is beginning of GI (CP) 



% construct the channel impulse response matrix 
for i=l:M 
for j=l:N 

if i-j+l<=0 

break; 
else 

H(i, j)=ch(i-j + l) ; 

end 

end 

end 
n=l; 

for i=M+l:M+N-l 
n=n+l ; 
k=0; 

for j=n:N 

k=k+l; 

H(i f j)=ch(M-k+l); a a 

end FIG. 4A 

end 



for i=l:mu 
\ mv(i)=l; 
\ end 

\ C= conv (mv, abs (ch) ) ; 
\ [Y, I] = max (abs (C) ) ; 



* 



for i=l:mu 
mv (i) =1; 
' end 

C= conv (mv, abs (ch) ) ; 
[Y,I] = max (abs (C) ) ; 
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initialize the SIRF filter to an arbitrary value 
:=[ -0.0593 
0.01047 
-0.062386 
0. 02418 
-0.065114 
0.030031 
-0. 039083 
0. 13789 
-0. 10266 
-0.014681 
-0. 11224 
0.0546 
-0. 12642 
0. 18608 
-0.020895 
0.38407 
-0.30117 
-0.37885 
0. 42326] ; 



/^g = f ft (c, length) ; 
S f = abs (g) ; 
7 lamda = 1 ; 
[ tol = 50; 
V, cold = sum (abs (c) ) ; 



FIG. 4B 



%begin iteration 

for ii=l:20 % 20 is the number of iterations 
% projection on the the set C2 
for m = wp: length/2 

if ( f (m) > (1+beta) ) 
gg=g (m) ; 

g (m) = (1 + beta) /f (m) ^complex (real (g (m) ), imag (g (m) )) ; 
g (m) = gg+lamda* (g (m) -gg) ; 
end 

if ( f (m) < (1-beta) ) 
gg=g (m) ; 

g(m) = (1 - beta) /f (m) ^complex (real (g (m) ), imag (g (m) ) ) 
g(m) = gg+lamda* (g (m) -gg) ; 

end 

end 

for m = l:ws 

if (f (m) > delta) 
gg=g(m) ; 

g (m) = delta/f (m) ^complex (real (g (m) ) , imag (g (m) ) ) ; 
g (m) = gg+lamda* (g (m) -gg) ; 

end 

end 

for m =■ 2 : length/2 

g (m+length/2) = conj (g (length/2-m+2 ) ) ; 

end 

cr = real (if ft ( (g) , length) ) ; ^transform to time domain 
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% time domain projection on C4 set 
for i = 1:N 

c(i) = cr(i); %project into the set C4 

end 

% time domain projection on C5 set 
for n=I-mu-10 : I-mu 
norm = 0; 
prod = 0; 
%prodl=0; 
for i=l:N 

norm = norm + H(n,i) A 2; 
prod = prod + c(i)*H(n,i); 
%prodl = prodl + c (i) *ch (n-i+1) ; 

end 

if (prod > tol) 
for nn=l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 

end 

end 

if (prod < -tol) 
for nn=l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (-tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 



end 

for n=I:I+100 
norm = 0 ; 
prod = 0; 
for i=l:N 

norm = norm + H(n,i) ^ 2; 

prod = prod + c(i)*H(n / i); 

end 

if (prod > tol) 
for nn=l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamda* (c (nn) -cc) ; 



end 



end 



end 

end 

if (prod < -tol) 
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for nn-l:N 

cc = c (nn) ; 

c(nn) = c(nn) + ( (-tol-prod) /norm) * H(n,nn); 
c(nn) = cc+lamcia* (c (nn) -cc) ; 

end 

end 

end 

ss = abs (cold-sum (abs (c) ) ) 
cold = sum(abs(c)); 

g = f ft (c, length) ; %transform to frequency domain 
f = abs (g) ; 

end 

for i = 1:N 

sirf(i) = c(i); %sirf will hold the SIRF filter coefficients 
end 
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